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ABSTRACT 

UV radiation from early astrophysical sources could have a large impact on subsequent star 
formation in nearby protogalaxies, and in general on the progress of cosmological reioniza- 
tion. Theoretical arguments based on the absence of metals in the early Universe suggest that 
the first stars were likely massive, bright, yet short-lived, with lifetimes of a few million years. 
Here we study the radiative feedback arising from such stars using hydrodynamical simu- 
lations with transient UV backgrounds (UVBs) and persistent Lyman-Werner backgrounds 
(LWBs) of varying intensity. We extend our prior work in Mesinger et al. (2006), by studying 
a more typical region whose proto-galaxies form at lower redshifts, z ~ 13-20, in the epoch 
likely preceding the bulk of reionization. We confirm our previous results that feedback in 
the relic HII regions resulting from such transient radiation, is itself transient. Feedback ef- 
fects dwindle away after ~ 30% of the Hubble time, and the same critical specific intensity 
of Juv ~ 0.1 x 10~ 21 ergs s _1 cm~ 2 Hz -1 sr _1 separates positive and negative feedback 
regimes. This suggests that overall feedback is fairly insensitive to the large-scale environ- 
ment, overdensity, and redshift-dependent halo parameters, and can accurately be modeled in 
this regime with just the intensity of the impinging UVB. Additionally, we discover a second 
episode of eventual positive feedback in halos which have not yet collapsed when their pro- 
genitor regions were exposed to the transient UVB. When exposed to the transient UVB, this 
gas suffers relatively little density depletion but a significant enhancement of the molecular 
hydrogen abundance, thus resulting in net positive feedback. This eventual positive feedback 
appears in all runs, regardless of the strength of the UVB. However, this feedback regime is 
very sensitive to the presence of Lyman-Werner radiation, and notable effects disappear under 
fairly modest background intensities of Jlw > 10 -3 x 10 _21 ergs s" 1 cm~ 2 Hz -1 sr" 1 , as- 
suming the region is optically thin for LW photons. Nevertheless, when exposed to the same 
LWB, halos inside relic HII regions always have a higher H2 abundance and shorter cooling 
times than halos outside relic HII regions, allowing gas to cool faster once it finally begins to 
collapse onto the halo. We conclude that UV radiative feedback in relic HII regions, although 
a complicated process, seems unlikely to have a major impact on the progress of cosmolog- 
ical reionization, provided that present estimates of the lifetime and luminosity of a Pop III 
star are accurate. More likely is that the build-up of the LWB ultimately governs the feedback 
strength until a persistent UV background can be established. 
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1 INTRODUCTION 

Semi-analytic models and numerical simulations predict that the 
first astrophysical objects formed at redshifts z > 20. They form 
out of metal-free gas with inefficient fragmentation, thus numeri- 
cal simulations predict that these first objects could be very mas- 
sive (M« IOOMq), so-called 'Population III' (Pop I II) stars 
dAbel et aljEool ; iBromm et alj|2002l ; [Yoshida et alj|2008h . These 
stars would be short-lived, with lifetimes of ~ 3 Myr, but could 
produce ten times more ionizing photons per baryon than regular 
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Population II stars JSchaereil2 002. 2003). Also, such star formation 
could have been somewhat synchronized, given that these early ob- 
jects formed in highly biased and clustered environments. Hence, 
Pop III stars could have a significant impact on subsequent genera- 
tions of objects. 

Indeed Pop III stars can strongly affect the progress of 
cosmological reioniza tion (jHaiman & Holder 2003; Cen 2 003bl l3; 
IWvithe & Loebi r2003). although their contribution is virtually un- 
known from first principles and is generally applied to simulations 
with a "tuning knob" approach. Thus ionizing photons from the 
first stars and la ck thereof have b een used to explain bot h early 
(e.g.|Cerj2003a) and fairly late (e.g. lHaiman & Brvanl20 06) reion- 
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ization, depending on which scenario was favored when the works 
were published. Therefore, it is quite important to gain physical 
insight into how the radiation of the first stars impacted their sur- 
roundings. 

Radiative feedback can be either positive or negative, in that it 
can enhance or suppress subsequent star-formation. Positive feed- 
back can result wh en the enhanced free -electron fraction from ion- 
izing photons (e.g. I Oh & Haiman 2002) or hydrodynamical shocks 
(Shapiro & Kang 1 19871) catalyzes the formation of molecular hy- 
drogen (H2). If the ionizing background "turns-off ' (resulting in 
so-called "relic" HII regions), H2 cooling can provide the domi- 
nant cooling channel at high-densities and low temperatures. Con- 
versely, negative feedback can result from heating by io nizing radi- 
ation which can photo-evaporate gas in low-mass halos ( Ef stathiou 
19921: iBarkana & Loeblll999t iGnedid bOQOt IShapiro et alj|2004l ; 
Diikst ra et al.l |2004 iMesinger & Diikstral 120081) . Also, an active 
background of Lyman-Werner (LW) radiation (with photon en- 
ergies in the 11.18-13.6eV range) can diss ociate H2, thus de- 
creas i ng the gas's cooling capabilities (e.g., lHaiman et al. J 1997 , 



200rilCiardi et alj2000l:|Machacek et alj200lMwise & Abell2007 : 
O'Shea & Normadl2008l) ~ 

Indeed, simulations find radiative feedback to be nuanced 
at very hi gh redshifts. Positiv e feedback dominates in flash ion- 
ized gas dO'Shea et alj l2005h . or ga s exposed to a weak tran- 
sient ultraviolet background (UVB) (Mesinger et al. 2006, here- 
after MBH0 6) such as might b e present close to the edges of 
HII regions dRicotti et alj|2002l : iKuhlen & Madad 120051) . On the 
other hand, negative feedback in relic HI I regions can occur in re- 
gions closer to the Pop III star (MBH06;|Susa & Umemurdl2006l : 
lAhn & Shapirdl2007l ; lYoshida et alj|2007l) . Furthermore, radiative 
trans fer effects can al s o impact the strength an d sign of feedback 
(e.g. Illiev et afll2005t ISusa & Umemurdhood : lAbel et alj I2OO7I ; 
Whal en et al.ll2008l : IWise & Abelll2008bl) ! Because of variations in 
halo mass, collapse redshift, central baryon density, stellar spec- 
trum, and distance to the ionizing source, it is difficult to accurately 
deal with the large parameter space governing photoevaporation. 

In MBH06, we attempted to tackle this issue with a statistical 
approach. By giving up on simulating "realistic" (as much as this is 
possible given present uncertainties) relic HII regions left behind by 
deceased Pop III stars, we were able to explore larger swaths of pa- 
rameter space and have a fairly large sample of second-generation 
objects. We did this by applying various uniform ultraviolet (UV) 
and LW backgrounds on our entire simulation box, and statistically 
comparing the resulting evolution against the fiducial run without 
radiation. Although this approach clearly does not realistically treat 
the photoionization around a particular source, it might provide a 
good statistical description of a ensemble of reionization studies. 
Clearly, at the onset of early reionization not all halos lie within 
HII regions, as they do in our models. Our primary purpose is not to 
simulate early reionization or the photoevaporation of well-formed 
Pop III halos in detail (which would require radiative transfer) but 
to evaluate the impact of relic HII regions on the formation of halos 
at later times. Our uniformly illuminated boxes provide an ensem- 
ble of evaporated halos at high redshift that statistically sample the 
effects of early photoionization on the assembly of halos and for- 
mation of cold dense gas at lower redshifts. 

However, since our analysis in MBH06 focused on a highly 
biased region which went non-linear at a high redshift, we could 
not follow the evolution to redshifts much lower than z ~ 20. 
This meant that we could only extrapol ate the observed trend s sug- 
gesting that the feedback was transient. IOh&Haimai] (b003) claim 
that, even in the presence of a weak LWB, excess entropy will even- 



tually suppress star formation in protogalaxies forming inside relic 
HII regions. In MBH06, we were unable to conclusively verify or 
refute this claim. 

In this work, we extend the study by MBH06 by analyzing 
radiative feedback in a less-biased region, whose typical proto- 
galaxies form at lower redshifts, z ~ 13-20. Given that the 5- 
yr WMAP polarization data place a constraint on the reionization 
redshift (assuming in stantaneous reionization) of z — 11.0 ± 1.4 
jDunklev et al.1 120081) . such cosmological regions are likely to 
host the majority of feedback effects from Pop III stars, before 
a persistent UVB a nd/or a strong LWB become entrenched (e.g. 
Haim anet al.ll 19961) . We apply the same statistical approach as in 
MBH06, and study various combinations of transient UVBs and 
persistent LWBs. As such, we attempt to provide a framework 
for analytically incorporating such radiative feedback into semi- 
analytic, and large-scale numerical and semi-numerical studies. 

In fj2]we enumerate and describe the simulations which are 
used in this work. In Sj3] we study the initial radiative feedback 
caused by a transient UVB, adding a persistent LWB in £13. 1 1 Then 
in Sj4] we describe a related mechanism that results in eventual pos- 
itive feedback, regardless of the strength of the UVB. In Sj5] we 
discuss the impact of neglecting self-shielding. Finally, we offer 
our conclusions in ^6] 

Throughout this paper, we adopt the background cosmologi- 
cal parameters (ft A , Qm, fif>, n, <r 8 , #0) = (0.7, 0.3, 0.047, 1, 0.92, 
70 km s -1 Mpc -1 ), consistent with the measurements of the power 
spectrum of CMB tempera ture anisotropics b y the first year of data 
from the WMAP satellite dSpergel et al]|2003l) . Alt hough the 5-yr 
data prefers a slightly low er value of us = 0.82 JDunklev et al.1 
bOOSt [Komatsu et ali boOS") which somewhat delays structure for- 
mation, we keep the cosmology the same as in MBH06, to facilitate 
direct comparison. Unless stated otherwise, we quote all quantities 
in comoving units. 



2 SIMULATIONS 

We use the Eulerian adaptive mesh refinement (AMR) code 
Enzo, which is describ ed in greater detail elsewhere ( Brvad [T999l : 
iNorman & Brvadll999h . Our simulation volume is 1 (ft 1 Mpc) s , 
initialized at Zj n j t = 99 with density perturbations drawn from the 
lEisenstein & Hd d 19991) power spectrum. Our root grid is 128 3 . We 
have two additional static levels of refinement inside a central 0.25 
ft -1 Mpc cube. In addition, grid cells inside the central region are 
allowed to dynamically refine so that the Jeans length is resolved 
by at least 4 grid zones and no grid cell contains more than 4 times 
the initial gas mass element. Each additional grid level refines the 
mesh length of the parent grid cell by a factor of 2. We allow for 
a maximum of 10 levels of refinement inside the refined central 
region, granting us a spatial resolution of 7.63 ft -1 pc. This co- 
moving resolution translates to 0.36 ft -1 proper pc at z = 20. 

As stated above, our simulation runs were set up to facilitate 
comparison with MBH06. Users interested in de tails of the sim- 
ulation s are encouraged to consult MBH06 and Machacek et al. 

There is one notable difference between our runs here and 
those in MBH06: our central refined region for these runs is not 
centered on the highest density region of the box. Instead is it cho- 
sen to be more typical of the regions expected to host the bulk of 
halos which form prior to any significant cosmological reioniza- 
tion. Specifically, the density inside our central (0.25ft -1 Mpc) 3 
corresponds to a 0.75 a mass fluctuation of an equivalent spherical 
volume (in MBH06, we studied a 2.4 a region). The less overdense 
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region allows us to extend the analysis of MBH06 to lower red- 
shifts, and test the robustness of the conclusions of MBH06 on a 
less-biased and more typical region. 

As shown in Table Q] we have performed four different runs 
without a LW background, distinguished by the duration or ampli- 
tude of the assumed UVB, and eight additional runs that include 
an additional constant LW background. Again, these runs are anal- 
ogous runs to those in MBH06. For the UV radiation we assume 
an isotropic background flux with a T = 2 x 10 4 K blackbody 
spectral shape, normalized at the hydrogen ionization frequency, 
hvH = 13.6 eV. This spectrum is softer than the ~ 8 x 10 4 K 
blackbody sp ectrum predicted to be produced by Pop III stars (e.g. 
ISchaererl2 002). which leads us to find a lower reionization temper- 
ature (by about a factor of two) than we would have found with a 
harder spectrum. This results in outflow velocities (and therefore 
halo evacuation times) which are also too low by this factor. 

The values of J^x are shown in Table Q] in units of 
10~ 21 ergs s _1 cm -2 Hz -1 sr _1 . The NoUVB run contains no 
UV radiation, and serves mainly as a reference run. The Heat0.08 
and Heat0.8 runs include a UVB with J^i = 0.08, 0.8, respec- 
tively. The value of J^ v = 0.08 was chosen to correspond to the 
mean UV flux expect ed inside a typical H II region surrounding a 
primordial star (e.g., lAlvarez et al] |2006). As we do not include 
dynamically expanding HII regions in our code, the Heat0.8 and 
Flash (in which the UVB is turned on for a negligibly short dura- 
tion) runs can be viewed as extremes, corresponding to conditions 
close to the center and close to the edge of the HII region, respec- 
tively. More generally, studying a range of values of Juv is useful, 
since the radiative properties of high-redshift sources are uncer- 
tain. In the latter two runs, the UVB is turned on at zuvB,on = 25 
and turned off at zuvb.oS = 24.64. This redshift range corre- 
sponds to a typical theoretical stellar lifetime, ~ 3 Myr, of a 
~ IOOA/q primordial (Pop-Ill) star jSchaereJl2002h . The flash 
ionization run, Flash, instantaneously sets the gas temperature to 
T=15000 K and the hydrogen neutral fraction to xm = 10~ 3 
throughout the simulation volume, but involves no heating there- 
after. Again, our adopted temperature is probably too low by a fac- 
tor of two compared to gas heated by Pop III stars. Finally, the eight 
runs in the bottom half of the Table repeat pairs of the NoUVB and 
Heat0.8 runs with four different constant LW backgrounds (J|i = 
10~ 4 , 0.001, 0.01, and 0.1, normalized at 12.87eV in units of 
10 _21 ergs s _1 cm -2 Hz -1 sr _1 , and assumed to be frequency- 
independent within the narrow LW band). 

We obtain the dark matter (DM) and baryonic mass of halos, 
in the same manner as MBH06. Similarly, we use the fraction of 
total gas within the virial radius which is cold and dense (CD), / c d 
as a proxy for the gas reservoir available to form stars. By cold, 
we mean gas whose temperature is < 0.5T v i r , where T v i r is the 
halo's virial temperature (for how viri al temperatures are a ssoci- 
ated with halos in the simulation, see lMachaceketal1 l200lh. By 
dense, we mean gas whose density is > 10 19 Mq Mpc~ 3 ss 330 
cm" , roughly corresponding to the density at which the baryons 
become important to the gravitational potential in the co re, taken to 
be an immediate precursor to primordial star formation dAbel et al.l 
l2002h . 

Additionally, we discount halos which have been substantially 
contaminated by the large (low-resolution) DM particles originally 
outside of our refined region which enter it as it becomes nonlinear. 
Specifically, we remove from our analysis halos with an average 
DM particle mass greater than ~140% of the refined region's DM 
mass resolution, 747 Mq . However, since there are small stocastic 
variations in the amount of contamination between the runs, we first 



Table 1. Summary of Simulation Runs 

Run Name J^ x v z UV B,on ^uVB.off ^21 



Runs without a LW background 



NoUVB 





NA 


NA 





Flash 




25 


25 





Heat0.08 


0.08 


25 


24.62 





Heat0.8 


0.8 


25 


24.62 





Runs with a LW background 


LWl()- 4 _NoUVB 





NA 


NA 


10- 4 


LW10- 4 _Heat0.8 


0.8 


25 


24.62 


10~ 4 


LW10~ 3 .NoUVB 





NA 


NA 


10~ 3 


LW10~ 3 _Heat0.8 


0.8 


25 


24.62 


10- 3 


LW10" 2 .NoUVB 





NA 


NA 


0.01 


LW10- 2 _Heat0.8 


0.8 


25 


24.62 


0.01 


LW0.1 -NoUVB 





NA 


NA 


0.1 


LW0.1.Heat0.8 


0.8 


25 


24.62 


0.1 



apply this cut-off in the NoUVB run, and then remove the same 
halos from the analysis of the other runs. This way, we study the 
same halos in each run. 

Another possible source of contamination arises from closely 
separated halos. If some CD gas belonging to a halo is within an- 
other halo's virial radius (most likely in the process of merging), 
the other halo could undeservedly be flagged as containing CD gas 
as well. To counteract this, we set / c( j = for all halos without 
CD gas in their inner-most spherical shell. CD gas invariably first 
forms close to the halo's core, where the cooling times are the short- 
est; thus halos which contain CD gas at large radii but none in their 
core are most likely in the process of merging with a halo with gen- 
uine CD gas. Note that this algorithm is different than the one used 
in MBH06, where we discriminated based on mass and separation 
between halos. In this work such contamination was more preva- 
lent and so we devised what proved to be a more robust means of 
discrimination. 

Finally, we note that we d o not include HP coolin g in 
these simulations. As is shown in iMcGreer & BrvarJ (2008) and 
lYoshida et all fe007h . HD cooling is important in primordial halos 
only in the density range 10 2 - 10 B cm , where H2 can bring the 
temperature below about 200 K. Therefore it does not change the 
amount of CD gas, our primary diagnostic (although it would lead 
to relatively small changes in the predicted masses of the resulting 
stars). 



3 TRANSIENT FEEDBACK FOLLOWING A TRANSIENT 
UVB 

We first wish to confirm whether our conclusions in MBH06 also 
apply to this new sample of halos, which form at later times in a less 
overdense region of the universe. In MBH06, we found the follow- 
ing. The UVB ionizes and heats the gas in our simulation box^Gas 
inside halos which used to be at or close to hydrostatic equilibrium 
feels an increased pressure gradient from the rise in temperature, 
resulting in an outward-moving pressure shock. Once the UVB is 

1 It is important to note that we are operating in the optically thin limit; the 
accuracy and impact of that approximation is discussed in section|3] 
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Figure 1. Spherically averaged radial profiles of the same halo in the NoUVB (solid lines), Heat0.08 (dashed lines) and Heat0.8 (dotted lines) simulation runs. 
This halo was first able to form CD gas at z = 20, 21, 19 in the NoUVB, Heat0.08, and Heat0.8 runs, respectively. The top and bottom pairs of figures show 
snapshots at z = 24.62, 23, 21, 19 (left to right). The halo mass in the NoUVB run increases from M=4.2 X 10 5 A/ Q at z = 24.62, to M=1.4 X 1O 6 M at 
z = 19. All quantities are shown in proper (not comoving) units. The Upper panels show the hydrogen density, mass-weighted gas temperature, gas cooling 
time, and radial velocity (clockwise from upper left). The Bottom panels show mass fractions of HI, HII, H2, and the number fraction of e~ (clockwise from 
upper left). 



turned off, the gas quickly cools to T ~ 10 3 K through a combi- 
nation of atomic, molecular hydrogen and Compton cooling. The 
molecular hydrogen formation time is shorter than the recombina- 
tion time and a large amount of molecular hydrogen is produced, 
Hh s ~ few x!0~ 3 , irrespec tive of the density and temperature 
(see also lOh & H aiman 2002)0 The pressure shock starts to dissi- 
pate and gas expelled from the center of the halo starts falling back 
in (for a more complete description, readers are encouraged to see 
MBH06). 

Quantitatively, we concluded that the net impact of a transient 
UVB on star-formation could be understood simply in terms of the 
molecular hydrogen cooling time-scale: 



£h 2 



1.5k B T 



Ah 2 nmnH 2 
1.5fc B T 1 



A H , 
\10 3 K J 



n g XH 2 
2.5 



(1) 
(2) 



XH 2 



3 x IO- 3 



lcm ~ 



Myr, 



Here, fcg is the Boltzmann constant, T is the temperature, Ah 2 is 
the H2 cooling function, x e is the free electron number fraction, n g , 
rim, and jih 2 are the number densities of all baryons and electrons, 
neutral hydrogen, and H2, respectively, and the second equality is 
accurate shortly after ZuVB,off> since the recombination times at 
high redshifts and inside halos are very short. 



2 iGlover &~A bel 2008 have recently shown that including cooling due to 
excitations of H2 through collisions with protons and electrons can lead to 
a lower molecular hydrogen fraction by a factor of 2. This would result in a 
longer delay in star-formation than shown here. 



We found that the net delay in star-formation caused by the 
transient UVB is simply the ratio of the H2 cooling time in the 
run with the UVB, to the H2 cooling time in the run without UV 
heating. This ratio is calculated near the halo core, shortly after the 
UVB disappears (when the shock first starts to dissipate and the 
gas behind it starts infalling again). For example, the delay in star- 
formation of a halo in the Heat0.8 run with respect to the same halo 
in the NoUVB run is: 



/delay 



.HeatO.8 

j-NoUVB 
r H 2 



ZyNoUVB 
I yHoat0.8 



HcatO.I 



(3) 

We will now demonstrate that eq. ([3j is equally adept at quan- 
tifying feedback in our new simulations. First, we show the spher- 
ically averaged radial profiles of the same (typical) halo in the 
NoUVB (solid lines), Heat0.08 (dashed lines) and Heat0.8 (dot- 
ted lines) simulation runs in our Figure Q] This halo, which has a 
mass of 4 x 1O 5 M at z = 24.62, was first able to form CD gas 
at z = 20, 21, 19 in the NoUVB, Heat0.08, and Heat0.8 runs, re- 
spectively. The top and bottom pairs of figures are snapshots at z = 
24.62, 23, 21, and 19 (left to right). The upper panels show the hy- 
drogen density, mass-weighted gas temperature, gas cooling time, 
and radial velocity (clockwise from upper left). The bottom panels 
show mass fractions of HI, HII, H2, and the number fraction of e~ 
(clockwise from upper left). 

The profiles confirm the qualitative story line from the begin- 
ning of this section. The outward-moving pressure shock is evi- 
dent in both UVB runs from the density and velocity panels at 
z = zijVB.off, just before the UVB was turned off. This gas has 
already cooled to T = 10 3 K and started to collapse back onto the 
halo only a few Myr afterward at z = 23. Note that the outflow 
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and the resulting suppression of gas near the core is stronger in the 
run with a stronger UVB. Halos in the runs with a UVB also ex- 
perience a strong boost in the H2 fraction following 2uvB,off, with 
the increase being insensitive to the strength of the UVB. Forma- 
tion of CD gas in this halo is notably delayed until z = 19 in the 
Heat0.8 run, while being mildly expedited in the Heat0.08 run. Can 
this behavior be understood through eq. l[3}? 

Let us focus on the output we have shortly after 2uvB,off, at 
2 = 23. We can note that at this redshift, the pressure-driven shock 
has almost dissipated and reached zero velocity with gas recom- 
mencing infall behind the shock. Plugging in appropriate values 
into eq.[3]from the inner annulus of the halo, we obtain values of 
/delay ~ 0.5 and 1.5 for the Heat0.08 and Heat0.8 runs, in good 
agreement with the simulation outcomes. This behavior again sup- 
ports our conclusions in MBH06, where the same analysis held for 
a different, randomly-chosen halo, and implies that we understand 
the physical picture associated with feedback in relic HII regions. 

In order to quantify the suppression of CD gas by a transient 
UVB on our entire sample of halos, we define (as in MBH06) the 
cumulative, fractional suppression of the number of newly formed 
halos as 



l(*) 



lyNoUVB I 
cd 



(4) 



(z) „ jV NoUVB (2UVBon) 

where N? d VVB {z) and N^ ni (z) are the total number of halos 
with CD gas at redshift z in the NoUVB run and some given run 
i, respectively. This expression is well-defined for N^ d ° VVB (z) 

^ a 7-NoUVB / \. f AfNoUVB/ \_ w NoUVB/ \ „ 

> iv cd (zuvB.ou), tor l\ cd (z)=l\ cd (zuvB.onJ, we 
set 5jv, cd (z) = 0. Note that by definition, A^ oUVB (zuvb.oi) = 

^d m (2UVB,on). 

Similarly, we define the cumulative, fractional suppression of 
the CD gas mass as 



5M,cd(z) 



where M„, 



'(ZOVB,. 



M NoUVB(. 



J. f NoUVB 



(zUVB,o 



1 



(5) 



s (z) and M c r d m (,z) are the total mass of CD gas 
at redshift z in the NoUVB run and some given run i, respec- 
tively. The total CD gas mass is obtained by merely summing 
the CD gas masses for all of the halos in the simulation. As for 
equation (|4j, this expression is well-defined for M^ oUVB («) > 



M NoUVB 



(zUVB.on), 



M c N d ° UVB ( ZuvB ,on), and for M? d ° VVB (z) 

we setd M ,cd( z ) — °- 

Equations {4} and $5$ provide an estimate of how the CD gas 
has been affected by the presence of a UVB, following the turn-on 
redshift of the UVB, ZuvB,on (the values at zuvB,on are subtracted 
in order to provide a more sensitive measure of relative changes of 
CD gas). As defined above, 5jv,cd(z) = and 8m,ca{z) = if 
the UVB has no effect. If the effect of a UVB is positive, resulting 
in positive feedback, (5jv jC d(z) and SM.cd(z) would be positive. If 
the effect of the UVB is negative, <5jv,cd(z) and SM,cd(z) would be 
negative. 

In Figure [2] we plot the values of M* d nl (z) (top panel) and 
N* d nl (z) (bottom panel) in our four simulation runs with no LW 
background: mm = NoUVB (crosses), Flash (dashes), Heat0.08 
(triangles), and Heat0.8 (squares). We note that the number of ha- 
los with CD gas increases in all runs until z ~ 15, when our cen- 
tral refined regions starts becoming more non-linear. At these lower 
redshifts, halos become contaminated with large DM particles en- 
tering from outside the refined region and we remove them from 
our analysis sample (see Sf2J. 

The corresponding values of 8m, c a{z) and 5n, c a{z) are plot- 
ted in Figure[3]in the top and bottom panels, respectively. Although 
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Figure 2. Values of M^ ni (z) (top panel) and A^ 1 ™ (2) (bottom panel) as 
denned in equations {5) and @. The results are displayed for the NoUVB 
(crosses), Flash (dashes), Heat0.08 (triangles), and Heat0.8 (squares) sim- 
ulation runs. 

some of the notable fractional changes shown in Figure [3] might 
appear statistically insignificant due to the small number statistics 
inferred from Figure [2] it should be noted that these runs are not 
uncorrelated experiments. In other words, each of our runs in Table 
Q]is seeded with the same initial conditions, and so small relative 
changes compared to the NoUVB run are significant (i.e. the errors 
are not Poisson). Nevertheless, there are some minor differences 
between the runs in the amount of "contamination" by DM parti- 
cles outside of our refined region. 

In MBH06, we presented two main inferences regarding the 
effects of a transient UVB (aside from the delay due to cool- 
ing times discussed above): (1) a "critical" background value of 
J^x ~ 0.1 separated positive from negative feedback (stronger 
backgrounds result in negative feedback, while weaker back- 
grounds result in positive feedback; see the discussion surround- 
ing eq. [3] above); (2) the feedback effects, regardless of the sign, 
were transient and started to "disappear" by z ~ 18. However, we 
were unable to evolve our simulations below z < 18, since the cen- 
tral refined region became too non-linear and further progress was 
computationally prohibitive. Thus point (2) was an inference based 
on a figure analogous to Fig. [3] 

In this work we present results from a less biased region and 
are able to go to lower redshifts, z ~ 13.5. Thus we can test both 
of the conclusions above and see if they are sensitive to the large- 
scale halo environment and redshift. From figure[3] we see that both 
conclusions are valid for these simulations as well. In other words, 
our transient (corresponding to the expected lifetime of a Pop III 
star) UVB results in transient feedback which disappears after ~ 
30% of the Hubble time, with J^x ~ 0.1 separating net positive 
and negative feedback. Thus, these conclusions seem to be fairly 
insensitive to the large-scale overdensity of the simulated region, 
and by extension to the redshift-dependent properties of halos at 

ZUVB.on- 
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Figure 3. Values of &M,cd( z ) ( t0 P panel) and &N,cd{ z ) (bottom panel) 
as defined in equations {5} and @. The results are derived from Figure [2] 
and displayed for the Flash (dashes), Heat0.08 (triangles), and HeatO.8 
(squares) simulation runs. 

It is important to note a new result evident in Fig. [3] all runs 
with a transient UVB experience eventual positive feedback at z < 
15. This is true regardless of the strength of the UVB. We will 
examine this interesting result further in Sj4] 

3.1 Adding a Persistent LW Background 

Next, we investigate how the combination of transient UV heating 
and a persistent LW background impact star formation. This is a 
more physically-relevant scenario, since a LW background is likely 
estab lished before a lar ge fraction of the Universe has been reion- 
ized faaiman et all2000l) . 

Photons in the LW band dissociate H2 molecules; thus by 
definition they provide additional negative feedback, undercutting 
some of the enhancement in the H2 abundances in relic HII re- 
gions. Negative feedback from a LW background kicks in when the 
H2 dissociation timescale becomes shorter than the H2 formation 
timescale. Since the formation timescale is inversely proportional 
to gas density, whereas the dissociation timescale is independent of 
density, the density decrease caused by UV heating should make 
halos more susceptible to t he negative feedback of a LW back- 
ground ( Oh & Haiman 2003)0. We investigate these processes fur- 
ther below. 

Specifically, in Fig. [4] we plot values of Sm,c<i(z) (top panel) 
and 5jv,cd(z) (bottom panel) for our runs including a LW back- 
ground. The LW background turns on at z = 24.62 and remains 
on. At the simplest level, we can see the suppression of molecu- 
lar hydrogen (and hence cooling) due to the LW background by 

3 Note that the temperature increase associated with a UVB does not con- 
tribute much this combined negative feedback because the ionized gas is 
able to Compton cool down to the virial temperatures of these minihalos. 



Figure 4. Suppression of cold dense gas in halos in simulations runs that 
include a persistent LW background. The LW background had specific in- 
tensities of J LW = 10~ 4 , 10~ 3 , 0.01. or 0.1 (normalized at 12.87 eV, in 
units of 10 —21 ergs s -1 cm -2 Hz -1 sr -1 ). Each of these three LW back- 
grounds is applied to both our NoUVB and Heat0.8 runs at z = 24.62, and 
is subsequently left on. Values of 8]\A,cd{z) (top panel) and <5jv,cd(-z) (bot- 
tom panel) are shown, as defined in equations (5} and @. 

the more negative values of both 5m, cd(z) and 5jv,cd(z) at fixed 
redshift. 

Does this panel show evidence for the additional negative im- 
pact of the transient UVB discussed above? Without UV heating, 
a LW background with a specific intensity of = 0.01 is 

needed for notable negative feedback (see the blue empty trian- 
gles at z < 14 in the lower panel). This value is similar to the one 
found in MBH06, where we showed that by equating Hb-cooling 
and Hb-photodissociation one expects this critical LW background 
to scale as J™ oc n s /H2/TQwhen the UV heating is added, this 
critical value of J21 causing negative feedback decreases by a fac- 
tor of ~ 10 to J\f ~ 10" 3 (see the green solid circles in lower 
panel at z < 14). This decrease seems to confirm the above ar- 
guments. However, such an interpretation is too simplistic and the 
negative impact of the LWB can be ameliorated by positive feed- 
back inside the HII region at lower redshifts, as we shall see in ij4] 

From Fig. |4]one can also note that a value of J™ ~ 10 -3 - 
10 -2 separates feedback regimes dominated by a LW background 
from those dominated by our transient UVB. This can be seen by 
the fact that the amount of suppression differs between the NoUVB 
and HeatO.8 cases at low values of the LW background (i.e. for 
■^2i V < 10 -3 ), while at large values of J21, the amount of sup- 
pression is independent of the UVB. Near the threshold value of 
Jji , negative feedback transitions from being UV heating dom- 
inated (at z > 17, or < 100 Myr after zuvB,off) to being LW 

4 Since the mean density of gas collapsing into halos in these simulations is 
less than a factor of two smaller than those in MBH06, we roughly expect 
this critical value of Ji^ to also be less than a factor of two smaller (at 
least in the linear regime). 
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background dominated (z < 17, or > 100 Myr after zuvB.off)- 
As highlighted in MBH06, this "transition" behavior can be un- 
derstood as the combined result of two effects: the UV heating is 
turned off, and its impact is transient, while the crit ical LW back- 
ground scales roughly with the inverse of the density dHaiman et al.l 
l2000l:IOh & Haimanll2003h and hence a fixed LW background will 
have a larger impact at lower densities or decreasing redshifts. 



4 EVENTUAL POSITIVE FEEDBACK 

So far, we have not focused on the interesting new result hinted 
by the previous figures (e.g. see the z < 14.5 points in Fig. [3}. 
This figure shows evidence of delayed (approaching our lowest 
redshifts) positive feedback in all runs with a UVB. This result is 
somewhat surprising given that the Heat0.8 run shows strong nega- 
tive feedback initially after zuvb.oS • In fact, as was already pointed 
out when discussing Fig. [3] the initial feedback, whether positive or 
negative, fades away in all runs by z ~ 17. 

In Fig. [5] we show the total gas fractions (top panels) and CD 
gas fractions (bottom panels) from halos at z — 13.5. The positive 
feedback mentioned above is immediately evident from this figure, 
and shows explicitly the masses of the affected halos. In particu- 
lar, the NoUVB run at z = 13.5 has 3 fewer halos hosting CD 
gas than the Flash and Heat0.08 runs, and 2 fewer than even the 
Heat0.8 run (see also Fig- [2j. These are some of the smallest and 
youngest of such halos at these redshifts. We also verify that these 
halos formed in less biased and less overdense regions than the ha- 
los which didn't experience such positive feedback. 

What is the physical cause of such delayed positive feedback? 
Before we attempt to answer this question, let us briefly review the 
sequence of feedback-related events arising from a transient UVB: 

(i) When the UVB turns on, gas gets ionized and heated to tem- 
peratures > 10 4 K. 

(ii) The temperature increase sets-off an outward moving pres- 
sure shock in the cores of halos, where density profiles have already 
steepened as gas approached hydrostatic equilibrium. This pressure 
shock smooths out the gas distribution and leads to a decrease in gas 
density in the cores of halos. 

(iii) Once the UVB is turned off, the gas rapidly cools to T ~ 
10 3 K through a combination of atomic, molecular hydrogen and 
Compton cooling. This temperature approximately corresponds to 
the gas temperature at the virial radius of such a proto-galactic, 
molecularly-cooled halo, thus effectively neutralizing the impact 
of temperature change on feedback. 

(iv) A large amount of molecular hydrogen is produced, xh 2 ~ 
few x 10~ 3 , irrespective of the gas density and temperature. 

(v) The pressure-shock begins to dissipate and gas with a newly 
enhanced H2 abundance starts falling back onto the partially evac- 
uated halo. 

It is the balance between the density depletion in step (ii) 
(which is dependent on the UVB intensity) and the enhanced cool- 
ing capability of the gas from step (iv) (which is independent of 
the UVB intensity), which determines the transient feedback out- 
comes discussed in § [3] But what happens to gas which is not in 
the cores of already-formed halos with a steepened density profile? 
In the extreme case of homogeneous gas, the UVB-induced tem- 
perature increase in (i) does not translate into a pressure shock. 
Thus, there is no density depletion from step (ii). In the more real- 
istic case of a clumpy gas distribution, pressure gradients will still 



form and smooth out the gas distribution, but the resulting deple- 
tion will be less than in the cores of halos in hydrostatic equilib- 
rium. However, such lower-density gas outside of halos will still be 
able to cool to temperatures of ~ 10 3 K since the Compton cooling 
time is independent of density, and at these high redshifts is over a 
factor of two shorter than the recombination time at mean density 
[in fact, these lower density regions take longer to recombine and 
hence can Compton-cool for a longer period reaching lower tem- 
peratures dOh & Haimani r2002): see also equations (2) and (3) in 
MBH06]. Likewise, step (iv) mentioned above also applies to such 
lower density gas, with the transient background resulting in an en- 
hanced molecular hydrogen abundan ce. Thus, we expect gas which 
has not yet formed a self-similar (e.g. lNavarro et aljl996h halo pro- 
file to experience weaker negative feedback [i.e. step (ii)], and the 
same amount of positive feedback [i.e. step (iv)], when exposed to 
identical transient UVBs. 

This conclusion is perhaps a bit counterintuitive, since 
many cosmological radiative feedback studies find that feed- 
back is a strong function of halo mass, with more massive, 



back ( Ef stathiou 


1992; 


Barkana & Loeb 


1999; iGnedinl [2000 


Shapiro et al. 2004; Kitayama & Ikeuchi 2000; Diikstra et al. 2004 


Mesinger & Diikstra 2008). However, there 


are two important dis- 



tinctions here: we are discussing feedback resulting from a tran- 
sient UVB; and we are comparing halos which have not yet formed 
and set up a self-similar gas profile at zuvB.off to those which have 
already formed, rather than comparing within each of these two 
groups. 

We also see some evidence of these trends in Fig. [6] where we 
plot the ratio (in runs with a UVB over the fiducial NoUVB run) 
of the H2 fractions (red squares), the temperatures (blue triangles), 
and the average gas densities (black circles) within the central 15pc 
(proper) of halos. These quantities are shown shortly after zuvB.off, 
at z — 23. As already noted in ij3]and in (iv) above, the H2 frac- 
tion tends to increase by the same factor, regardless of the nature of 
the heating. As a result, the sign of the overall feedback (negative 
or positive) is determined primarily by whether or not the depres- 
sion of the gas density overcomes the positive effects of this H2 
enhancement. The density depletion does depend strongly on the 
type and amount of heating. This depletion, shown in black circles, 
is of course stronger for the Heat0.8 than the Heat0.08 run. Thus ha- 
los in the former initially experience strong negative feedback and 
in the latter negligible or weak positive feedback, as we've seen in 
Sj3]and as noted by the green crosses in the figure. Note however 
that the difference in the amount of gas depletion between these 
two heating runs decreases with decreasing halo mass as one be- 
gins being affected by Jeans/filtering scale effect^ Following the 
above argument, this trend should continue and density depletion 
should decrease in regions which are less evolved (i.e. the black 
points should approach unity as the halo mass is decreased below 
the values shown in the figure). It is these regions which become 
the progenitors of our halos experiencing eventual positive feed- 
back. Unfortunately, at z — 23, these regions are too unevolved to 
be identified by our halo finding algorithm and thus do not appear 
in Fig. [6] 

This eventual positive feedback should be quite sensitive to 
the LWB. Indeed, looking at the bottom panel of Fig. [4] we see 



5 Note that accurately modeling Jeans smo othing requires very high reso- 
lution in both dark matter and gas iNaoz et al] 2009). These small scales are 
shown in Fig.[6]merely to show qualitative trends. 
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Figure 5. Total gas fractions (top panels) and cold, dense gas fractions (bottom panels) as a function of total halo mass at redshift z = 13.5, the lowest 
redshift output of our simulations. The four panels correspond to the NoUVB (top left), Flash (top right), Heat0.08 (bottom left), and Heat0.8 (bottom right) 
simulation runs. Note that all runs with a UVB have more halos with CD gas than the NoUVB run. 



that this feedback effect is wiped-out if J™ > 10 3 . It is 
likely that such modest values of J^i are reached well b efore the 
majority of the Universe is ionized taaiman et aTl feoOCn. though 
the buildup of the LWB can be slowed by an increasing opti- 
cal depth to LW photons from the en hanced H2 abundances in 
relic HII re gions |joh nson et al. 2007 ), or in shells s urrounding 
HII r egions feicotti et a l. 2001 ; Kuhlen & Madaj2005l : though see 
iKitavama et al .1 120041) . Thus studying the true physical importance 
of this positive feedback mechanism would have to be done with 
self-consistent cosmological simulations of reionization. 

Now, let us look at specific examples by studying the evolu- 
tion of one such halo experiencing eventual positive feedback. In 
Figure [7] we plot the same quantities as in Fig. Q] but for a halo 



which experiences positive feedback in the Heat0.8 run. Note that, 
unlike the more evolved halos which are subjected to the transient 
feedback and have already formed at ztjvb.oS, the earliest snap- 
shot we have of this halo is at 2 = 21. Solid and dashed curves 
correspond to the NoUVB and Heat0.8 runs, respectively. 

It is important to note that this halo forms considerably af- 
ter zuvB,off, and thus the "density depletion" seen in the figure 
at z > 17 is of a different nature than discussed in [J3]and Fig.[T] 
Namely, this gas wasn't evacuated from the halo center with a pres- 
sure shock. Instead, these late forming halos start pulling in baryons 
when their virial temperature exceeds that of the gas temperature. 
The heating and subsequent Compton cooling following the tran- 
sient UVB sets the IGM gas temperature in these runs at T ~ 10 3 
K. Thus baryonic infall onto newly formed halos can commence 



© 0000 RAS, MNRAS 000, 000-000 



High-z Relic HII Regions and Radiative Feedback 9 



z=23 



z = 23 




0.01 



10 2 - 



10 



0.1 



0.01 





□ 

□ 

□□ □ 


x 

□ 


Heat0.8 : 






□ 

□ 


- 






n 




: a A a 

: «* A 




\* * 




: & o °S= ' o co 

: %*p ° o 


O o 

O x 

X 


o 

o 
o 

o 


* T 

o n s 



10 5 



10 7 



M/M s 



M/M s 



Figure 6. Ratio of the H2 fractions (red squares), the temperatures (blue triangles), and the average gas densities (black circles) within the central 15pc 
(proper) of halos. Ratios of these quantities in the Heat0.08 and Heat0.8 runs over the fiducial NoUVB run are shown in the left and right panels, respectively, 
at z = 23. The delay factors from eq. ^3) are also plotted as green crosses. 



when T vir > 10 3 K. Since the IGM is cooler in the NoUVB run, 
gas there gets a "head-start" in accreting onto DM halos. However, 
as is evident in the figure, the enhanced H2 cooling channel in the 
Heat0.8 run allows the gas to quickly catch-up and surpass its twin 
in the NoUVB ran. Thus by redshift 2 = 14, the halo has formed 
CD gas in the Heat0.8 run, but not in the NoUVB run. 

How does adding a LWB impact this positive feedback? From 
Fig. [4] we have already surmised that a background intensity of 
■^21^ > 10~ 3 , can offset these trends by decreasing the en- 
hanced H2 fraction in relic HII regions. In Fig. [8] we plot pro- 
files from the same halo in the LW10" 3 .NoUVB (solid curves) 
and LW10~ 3 JTeat0.8 (dashed curves'). By comparing the solid 
curves from Fig. [7] to the dashed curves from Fig. [8] one confirms 
that a LWB with J^i* ~ 10 -3 effectively neutralizes the posi- 
tive feedback from the transient UVB. Thus at 2 = 14 the halo 
core has a similar gas profile and cooling time in the NoUVB and 
LW10~ 3 _Heat0.8 runs. 

By comparing the two profiles at z = 14 in Fig. [8] we see 
that the transient UVB still stimulates positive feedback, but the 
effect is not nearly as dramatic. Specifically, in the halo core at 
z = 14 the density is enhanced and the cooling time is lowered 
(both by a factor of ~ 2) in the LW1CP 3 _Heat0.8 run compared 
to the LW10" 3 _NoUVB run. This level of positive feedback is 
much more modest compared with the case with no LWB, where 
the equivalent density enhancement was a factor of ~ 50 and the 
cooling time was lower by a factor of ~ 15 (see Fig.|7J. 

Does this trend continue as we increase the strength of the 
UVB? In Fig. |9] the solid and dashed curves correspond to the 
LW10" 2 JMoUVB and LW1(T 2 _Heat0.8 runs at z = 14, respec- 
tively; similarly in Fig. [TO] the solid and dashed curves correspond 
to the LW0.1 .NoUVB and LW0.1.Heat0.8 runs at z = 14, respec- 
tively. Here we see that the core density in the LW0.1_Heat0.8 run 
is lower than in the LW0.1_NoUVB run, but the cooling time is 
also lower by the same factor of ~ 2 seen in Fig. [8] sourced by 



the enhanced H2 fraction. In fact, from higher redshift output and 
the previous figures, we see that the density is still "catching up" 
as time progresses; namely, the dashed density curve keeps getting 
closer to the solid density curve. Thus we see that halos inside relic 
HII regions, when exposed to the same LWB, always have a higher 
H2 abundance and shorter cooling times than halos outside relic 
HII regions, allowing gas to cool faster once it finally begins to 
collapse onto the halo. 

It is important to note that in this regime where a strong LWB 
dominates the overall feedback, star formation is likely to be sub- 
stantially delayed. This delay can be long enough to allow the halo 
to be photoevaporated during reionization, before it had a chance to 
host additional stars. Thus if one is only concerned with the fate of 
those minihalos contributing their ionizing photons to the progress 
of reionization, it is possible that the sole persistent feedback mech- 
anism is positive. 



5 THE IMPACT OF NEGLECTING RADIATIVE 
TRANSFER 

As we have noted a number of times in this paper, our simulations 
assume the optically thin limit, for both ionizing radiation and pho- 
todissociating LW photons. Both of these approximations clearly 
fail at points in our simulations, and so it is necessary to understand 
the nature and implications of this approximation. For clarity, we 
discuss ionizing photons and LW photons separately. 

We discuss the impact of radiation hydrodynamic (RHD) 
effects due to ionizing photons first. There is now a substan- 
tial body of literature o n the ionization of primordial gas ha- 
los dShamro et alj I200I: Iwhalen et alj |20o4; llliev et all | 2005|; 
I Alvarez et al. II2006I; ISusa & Umemurj|2006l ; Iwhalen et all [2008 ; 
lAbel et al.l2007l ; IWise & Abell2008bh and the broad outline of how 
halos are ionized is now well understood. We begin by noting that 
the optical depth of the halos which are ionized in our simulation 
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Figure 7. The same quantities as in Fig.[T] but for a halo which experiences positive feedback in the Heat0.8 run. Solid and dashed curves correspond to 
NoUVB and Heat().8 runs, respectively. 
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Figure 8. Same as Fig. [7] but with solid and dashed curves corresponding to the LW10 ' 5 _NoUVB and LW10 3 _Heat0.8 runs, respectively. 



box at z = 25 (as shown in Fig. Q3, as measured from the virial 
radius to the center is - at the Lyman limit - well over 100. There- 
fore, at first blush, the optically thin limit seems like a poor approx- 
imation. However, this impression is not entirely accurate for two 
reasons. 

First, we are studying the later formation of halos in relic HII 
regions, and much of the gas that eventually falls into a halo at 
late times is not in the center of h alos at reionization. This is the 
question that|5h]& Haimanj d2003h asked: does the excess entropy 
generated in low-density gas in relic HII regions prevent it from 
forming stars at later times. This low density gas has much lower 



optical depth and so is very likely to be ionized by a nearby Pop III 
star. 

Second, even for the gas in the centers of our halos at zuvB.on, 
the optically thin approximation is not as bad as the large optical 
depth at the Lyman limit indicates. We can see this more clearly 
by comparing our results to RHD calculations. For example, take 
the (typical) halo shown in Fig. [TJ The central density when ion- 
ization begins (at z = 25) is just slightly over 1 cm -3 . This 
value is much below the central density (10 4 cm -3 ) assumed by 
ISusa & Umemural d2004l) . and is not that different from their crit- 
ical density for self shielding (depending on the assumed distance 
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Figure 9. Same as Fig. [7] at z = 14, but with solid and dashed curves cor- 
responding to the LW10 _2 JSfoUVB and LW10 _2 _Heat0.8 runs, respec- 
tively. 



of the source). A better comparison is with the "023" series of sim- 
ulations run in IWhalen et al. (2008), which have almost identical 
central densities. As is shown in Fig. 17 of that paper, this halo is 
disrupted and dissociated for all assumed distances of the first star. 
Comparing Fig. 9 and 10 of that paper to our Fig. Q] shows some 
obvious differences due to radiative transfer effects, but the basic 
evolution is not all that different. In particular, after a few Myr, the 
density at the halo center has decreased by about a factor of ~ 5 af- 
ter ~ 3 Myr, while the molecular hydrogen fraction has increased 
to a few times 10 -3 . These numbers are very similar to what we 
find, and it is these basic facts which, we argue, cause the delay in 
future star formation, as discussed earlier. The reason for the sim- 
ilarity is that, as it sweeps over the halo, the ionization front can 
quickly ionize hydrogen until it transitions to a D-type front and 
stalls. However, this happens only ~ 15 pc from the core, at densi- 
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Figure 10. Same as Fig. [7] at z = 14, but with solid and dashed curves cor- 
responding to the LWO. 1 _NoUVB and LWO. 1 _HeatO. 8 runs, respectively. 

ties slightly below 1 cm , and the shock generated by the D-type 
front quickly disperses the core. 

Of course, we do not claim to reproduce the radiation hy- 
drodynamics in detail. For example, we do not repr oduce the thin 
self-sh ielded filament at the lower right of Fig. 10 of IWhalen et al.l 
(2008), nor do we follow the detailed shock-compression or 
"rocket" offset of the core. However, by choosing to apply the ion- 
izing flux at a redshift when the halos have low masses and cen- 
tral densities, we minimize the impact of the radiation transfer. To 
verify that this is actually true in our simulations, we examined 
the central densities and masses of our halos at z — 25; we find 
that all halos except for one have masses below 10 6 Mq and cen- 
tral densities below ~ 2 cm" 3 (one halo has Al > 1O 6 M0 and 
n c — 7 x 10 3 cm -3 ). Our halo abundance is consistent with pre- 
dictions from the conditional Press-Schechter formalism (MBH06). 
We chose this approach because we are interested in investigating 
the impact of relic HII regions on the formation of structure which 
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forms significantly later, rather than studying the impact of a Pop 
III star on a nearby, highly collapsed halo. In particular, as noted in 
th e introduction, our m otivation is to follow up on the suggestion 
of lOh & Haimanl f2003h that relic HII regions will imprint an ex- 
cess entropy on the low-density gas at early times, which will then 
suppress later star formation. Of course, we stress that it is still im- 
portant to examine our results with a more realistic RHD treatment. 

The results discussed in the section on positive feedback (sec- 
tion are even more robust because in that case the halos form 
out of gas that was barely overdense when the ionization occurred. 
As can be seen in Fig. [7] even at z — 21 (long after the transient 
ionization), the central density of the halo (in the noUVB case) is 
only 0.2 cm -3 . 

Another useful comparison of our results can be made with 
IWise & Abell |2008a), who used a full cosmological simulation of 
Pop III star formation including radiative transfer (albeit in a re- 
gion considerably smaller than considered th e re). W hile a direct 
comparison is difficult because IWise & Abeil (2008a) do not ex- 
plicitly measure the suppression due to photo-ionization, a number 
of conclusions can be drawn. First, it is clear from Fig. 7 of that 
paper that they also see the evaporation of gas from low-mass halos 
and filaments, and the resulting density field in the two simulations 
looks remarkably similar. In addition, the time difference between 
the first Pop III object and the formation of the second is about 
80 Myr, consistent with the time-delay derived in this paper. On 
the other hand, one important point that comes out of that paper 
is that, because Pop III star formation tends to be clustered, relic 
HII regions do not sit isolated for longer periods of time and so a 
region of space may undergo multiple ionization episodes, some- 
thing which is not considered in this work. 

Next, we turn to LW photons. First, we remind the reader that 
we are discussing a persistent LW background, and that although 
the LW flux is produced by the same class of objects that generate 
the transient UV flux, it is not produc ed by the same phy sical ob- 
jects. As discussed in more detail in lHaiman et alj dl996h . the LW 
background is accumulated from a large number of distant sources. 
The effective optical depth in the LW lines in t he background IGM , 
over a Hubble dis tance, is at most r ~ 1 — 2 (Haima n et alj|2000l ; 
iRicotti et alj|200l1) . so the Olbers' integral for the LWB is domi- 
nated by the large number of sources at a good fraction of the Hub- 
ble distance. In comparison, the ionizing photons come exclusively 
from the much fewer source(s) nearby, within the local HII region. 

LW photons can be self-shielded by a sufficient column of 
H2. Here again, there are two regimes. One is the large quantities 
of low-density gas in the relic HII regions that form substantial 
amounts of mole cular hydrogen immediat ely after the ionizing flux 
is turned off (e.g. lKuhl en & Madau 2005). Although H2 formed in 
this manner may slow the build-up of a persistent LW background, 
it does not - for a given LW background - pro vide a large amount of 
shielding. For example. I Johnson et al.1 12007) find typical values of 
LW optical depths of order unity for a single relic HII region. This 
means that the LW flux inside the relic HII region can be reduced 
by a factor of a few, but it also means that it doesn't take long for 
the persistent LW background to "burn off the H2 produced in 
relic HII regions. 

The second regime occurs after the ionized gas falls back into 
the collapsed halos, where H2 is in equilibrium between formation 
and photo-destruction. In this case, when a sufficient column den- 
sity of H2 is formed, it can effectively attenuate the LW flux. Based 
on our simulations, we find that when CD gas forms, the LW optical 
depth can be of order unity or larger, indicating that self-shielding 
can play a role during the recollapse (although probably not during 



the photoionization process itself). Therefore, we caution readers 
that the critical LW values we give above may need to be increased 
when self-shielding is properly accounted for. 

Finally, we note that in simulations with full radiative transfer, 
the photoionizing flux can also vary across a given halo (relic H 
II regions have radii of 2-3 kpc, while the Lagrangian radius of 
a typical halo in our simulations is of order 0.5 kpc). This also 
means that, while most forming halos lie either completely inside 
or outside a relic H II region, some proto-halos could be partially 
ionized. 



6 CONCLUSIONS 

UV radiation from early astrophysical sources could have a large 
impact on subsequent star formation in nearby protogalaxies, and 
in general on the progress of cosmological reionization. Theoretical 
arguments based on the absence of metals in the early Universe sug- 
gest that the first stars were likely massive, bright, yet short-lived, 
with lifetimes of a few million years. Here we study the impact of 
the transient radiation arising from such stars on early protogalax- 
ies. We apply the same statistical approach as in MBH06, study- 
ing various combinations of transient UVBs and persistent LWBs, 
using the hydrodynamical simulation code, Enzo. We also study 
a more typical and relevant region whose proto-galaxies form at 
lower redshifts, z ~ 13-20. This allows us to trace feedback ef- 
fects on longer time-scales and lower redshifts. 

We confirm the results of MBH06 that feedback in the relic 
HII regions resulting from such transient radiation, is itself tran- 
sient. Feedback effects dwindle away after ~ 30% of the Hub- 
ble time, and the same critical specific intensity of Juv ~ 0.1 x 
10~ 21 ergs s~ x cm -2 Hz -1 sr _1 separates positive and negative 
feedback regimes. A weaker UVB stimulates subsequent star for- 
mation inside the relic HII regions, by enhancing the molecular hy- 
drogen abundance. A stronger UVB delays star-formation by re- 
ducing the gas density at the centers of collapsing halos. The fact 
that we have confirmed the results in MBH06 (which were mostly 
inferences obtained by extrapolating from higher redshifts since we 
were unable to continue those simulations to low enough redshifts) 
suggests that overall feedback is fairly insensitive to the large-scale 
environment, overdensity, and redshift-dependent halo parameters, 
and can accurately be modeled in this regime with just the intensity 
of the impinging UVB. 

We note that a value of J™ ~ 10~ 3 -10~ 2 separates feed- 
back regimes dominated by the LW background from those dom- 
inated by the transient UVB. As discussed in section [5] there is 
some uncertainty in these values due to our optically thin treat- 
ment. In particular, accounting for LW self-shielding may increase 
the values of jji quoted here. 

Additionally, we discover a second episode of eventual pos- 
itive feedback in halos which have not yet collapsed when their 
progenitor regions were exposed to the transient UVB. Gas in 
such late-forming objects did not have time to set-up a steep den- 
sity profile, which would translate into a strong pressure shock 
when photo-heated. Instead they formed out of material which 
was of more moderate densities at zuvB,on where the lasting im- 
pact of the transient UVB was positive: seeding the gas with an 
enhanced abundance of molecular hydrogen, which aids in its 
eventual cooling. However, this feedback regime is very sensi- 
tive to the presence of Lyman- Werner radiation, and disappears 
under fairly modest background intensities of Jlw > 10~ 3 x 
10 _21 ergs s _1 cm -2 Hz -1 sr _1 . Nevertheless, when exposed to 
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the same LWB, halos inside relic HII regions always have a higher 
H2 abundance and shorter cooling times than halos outside relic 
HII regions, allowing gas to cool faster once it finally begins to col- 
lapse onto the halo. Although it is difficult to accurately estimate 
the build-up of the LWB, it seems likely that such modest values 
were surpassed well before the bulk of reionization dHaiman et al.l 
2000). Thus this persistent eventual positive feedback might only 
be of academic interest, provided LW self-shielding is not efficient. 

We conclude that radiative feedback from a short-lived UVB 
seems unlikely to have a major impact on the progress of global 
cosmological reionization, provided that present estimates of the 
luminosities a nd lifetimes of Pop III stars are accurate (e.g. 
ISchaeredl2002h . Thus it is likely that the LWB plays the dominant 
role in regulating feedback in relic HII regions. 
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